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We develop a non-perturbative method to calculate the 
density of states (DOS) p(io) of the fluctuating gap model de- 
scribing the low-energy physics of electrons on a disordered 
Peierls chain. For a real order parameter field we calculate 
p(0) (i.e. the DOS at the Fermi energy) exactly as a func- 
tional of the disorder for a chain of finite length L. Averaging 
p(0) with respect to a Gaussian probability distribution of the 
fluctuating Peierls order parameter, we show that for L — ► oo 
the average (p(0)) diverges for any finite value of the correla- 
tion length above the Peierls transition. Pseudogap behavior 
emerges only if the Peierls order parameter is finite and suf- 
ficiently large. 
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At low temperatures many quasi one-dimensional con- 
ductors become unstable and develop long-range charge- 
density- wave order, i.e. undergo a Peierls-transition p|. 
Within a mean field picture a finite value Ao of the 
Peierls order parameter leads for frequencies |a>| < |A | 
to a gap in the electronic density of states (DOS) p(u) 
||]. To achieve a better understanding of the effect of 
fluctuations on the Peierls transition, Lee, Rice, and 
Anderson || introduced the so-called fluctuating gap 
model (FGM). In this model the fluctuating part A(x) = 
A(x) — Ao of the order parameter is approximated by a 
Gaussian stochastic process with covariance 



(A(x)A(x')) = K(x,x') = K e 



(1) 



Here (. . .) denotes averaging over the probability distri- 
bution of A (a;), Kq is a positive constant, and £ is the 
order parameter correlation length. We assume that the 
field A(x) is real, corresponding to a charge-density wave 
that is commensurate with the lattice. 

Twenty years ago Sadovskii pj found an apparently ex- 
act algorithm to calculate the average DOS of the FGM. 
His calculations showed that for temperatures above the 
Peierls transition, in a regime where £ is large but finite, 
the average DOS exhibits a substantial suppression in the 
vicinity of the Fermi energy, a so-called pseudogap. The 
algorithm constructed by Sadovskii has also been applied 
in a different context to explain the weak pseudogap be- 
havior in the underdoped cuprates (5|. However, recently 
it has been pointed out |],|7] that Sadovskii's algorithm 
contains a subtle flaw and hence does not produce the 
exact DOS of the FGM. It is therefore important to com- 
pare this algorithm with limiting cases where (p(u>)) can 
be calculated without any approximation. 

Besides the limit £ — > oo where Sadovskii's algorithm 



is indeed exact , there exists another non-trivial limit 
where the exact (p(ui)) is known: if in Eq.(Q) we let 
£ — > 0, Kq — * oo, with _Ko£ — > D = const, the right- 
hand side of Eq.(|j]) reduces to 2D5Jx — x 1 ). As shown by 
Ovchinnikov and Erikhman (OE) [pf, in the limit L — > oo 
(where L is the length of the chain) the exact (p(oj)) 
can then be obtained from the stationary solution of a 
Fokker-Planck equation |j) . For small ui and Ao = one 
finds H {p{ijj)) oc \oj\n 3 |cu 1 1 1 . Singularities of this type 
at the band center of a random Hamiltonian have been 
discovered by Dyson ficfl , and have recently also been 
found in one-dimensional spin-gap systems JlT,^2|. It is 
important to note that in the FGM the singu arity is a 
consequence of the charge conjugation symmetry of the 
underlying Dirac Hamiltonian, and is not related to con- 
crete probability properties of A(x) J|[l3|. In particular, 
the singularity is not an artefact of the exactly solvable 
limit £ — > considered by OE ||. It is therefore reason- 
able to expect that for any £ < oo the average DOS of the 
FGM exhibits a singularity at u> = 0. This general argu- 
ment is in disagreement with Ref . Q , where for large but 
finite £ a pseudogap (and hence no singularity) has been 
obtained. In this work we shall resolve this contradiction 
by calculating the average DOS at the Fermi energy (i.e. 
{p{ijJ = 0))) exactly for arbitrary £. 

The local DOS p(x,uj) of the FGM for a given realiza- 
tion of the disorder can be written as 

p(x, u) = -n^ 1 lmTT[a 3 G(x, x, uj + i0)] , (2) 

where the 2x2 matrix Green's function Q satisfies 

[id x + u>a 3 — iA(x)a2) Q{x, x', to) — S(x — x')o~o . (3) 

Here o~i are the usual Pauli matrices and ao is the 2x2 
unit matrix. Note that in Eq.(||) we have factored out 
a Pauli matrix 173, so that the differential operator id x 
in Eq.(|3|) is proportional to the unit matrix. To solve 
Eq.(^|), we try the ansatz (suppressing for simplicity the 
frequency label) 



Q{x,x') = U{x)g l {x,x')U- 1 {x') , 



(4) 



where U(x) is an invertible 2x2 matrix. Eq.(^|) resembles 
the transformation law of the comparator in non-Abelian 
gauge theories p4| . In fact, Eq.(Q) can be viewed as 
a gauge transformation which generalizes the Schwinger 
ansatz 15 to the non-Abelian case. It is easy to show 
that the solution of Eq. (||) can indeed be written in the 
form (|^) provided Q\ and U satisfy 

[id x + uja 3 ] Qi{x, x ) — 5{x - x')a a , (5) 
id x U(x) = lu[U(x)ct 3 - <t 3 U(x)} + iA{x)a 2 U(x) . (6) 
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Eq.(g) defines the Green's function of free fermions, and 
can be solved trivially via Fourier transformation. The 
difficult part of the calculation is the solution of the ma- 
trix equation (||). We parameterize U(x) as follows, 



U[x) 



(7) 



where <r± = h[o~i i i& 2], and the three functions &±(x), 
$3(0:) have to be chosen such that U(x) satisfies Eq.([3|). 
A parameterization similar to Eq.(^) has recently been 
used by Schopohl JfJ) to study the Eilenberger equations 
of superconductivity. We find that the ansatz (^) solves 
Eq.(g) if $±(x) and $3(2;) satisfy 



d x $ + = -2iuj$ + + A(x)[l 



2iwi>_ - A(x)[l - 2$+$_ 



<9x* 3 = -iA(x)$ 



(8b) 
(8c) 



Non-linear differential equations of the type (|Sa|j are 
called Riccati equations. The set of equations obtained 
by Schopohl llfll h as a similar structure but is not iden- 
tical with Eqs.( pa[ - |8c1) . Note that Eq.(jBa|) involves only 
<f> + . If we manage to obtain the solution $ + , Eqs.( [3b| , |Sc| ) 
become simple linear equations which can be solved ex- 
actly. From Eqs.(|],[|) and (Q) it is easy to see that the 
local DOS can be written as 



p(x,uj) = n 1 Rei?(a;, u + iO) , R = 1 



2$+$_ . (9) 



Thus, to calculate the average DOS we have to average 
the product over the probability distribution of 

the field A(x). In the limit where the right-hand-side 
of Eq.(Q) reduces to 2DS(x — x') we can use the fact 
that $+ and $_ satisfy first order differential equations 
to express the average ($+$_) in terms of the solution 
of two coupled one-dimensional Fokker-Planck equations 
Jl7| |. Thus, our method leads to an algorithm for ob- 
taining the exact (p(u>)) without using the node counting 
theorem [||. However, the Fokker-Planck equation ob- 
tained by OE H within the phase formalism || is easier 
to solve than our system of two coupled Fokker-Planck 
equations. Hence, for J-function correlated disorder our 
approach does not have any practical advantage. 

On the other hand, if the disorder is not 5-function 
correlated, probability distributions of physical quanti- 
ties do in general not satisfy Fokker-Planck equations, 
and it is not so easy to perform controlled calculations or 
even obtain exact results |T^] . We now show that for real 
A(x) the local DOS at the Fermi energy can be calcu- 
lated exactly. To derive this result, let us introduce the 
complex vector 



-V2(l - 
r = I 1 - 2<!> . <i> 

V2$_ 




(10) 



Note that by construction R 2 — 2Z + Z^ — 1 for all x, and 
that according to Eq.(^|) the second component of ip is 
related to the local DOS. Using Eqs.(pa|,pb|) we find 



d x if> = -H(x)i/3 , H(x) = 2iujJ 3 + 2A(x)J 1 , (11) 
where Ji are spin J = 1 operators in the representation 
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1 
J 3 = ( 
0-1 



/, =-L I 1 1 I (12, 
v2 \ 1 



Eg. (|1 1[) is a linear multiplicative stochastic differential 
equation |L8] . Formally this equation looks like the imag- 
inary time Schrodinger equation for a J = 1 quantum 
spin in a random magnetic field, with x playing the role 
of imaginary time. Although the operator H in Eq. (|ll| ) is 
not hermitian, we may perform an analytic continuation 
to imaginary frequencies (lu = iE) to obtain a hermitian 
spin Hamiltonian. We thus arrive at the remarkable re- 
sult that the average DOS of the FGM can be obtained 
from the average state-vector of a J = 1 spin in a random 
magnetic field. Our non-linear transformation (|l(]) is well 
known in the quantum theory of magnetism: with the for- 
mal identification R — ► J3, \/2Z± — > ^J±, \/2$- — > tf, 
and \/2Q+ — > b, Eq.(|Tc|) is precisely the Dyson-Maleev 
transformation |19| , which expresses the spin operators 
in terms of boson operators b,b^ . 

The solution of Eq.([ll"|) with initial condition tp(0) = 
■00 is ip(x) = S(x)tpo, where the S-matrix is 



S(x) — Texp 



dx'H(x') 



(13) 



Here T is the usual time ordering operator. The proper 
choice of boundary conditions requires some care. For 
simplicity, let us assume that A(x) is non-zero only in 
a finite interval < x < L. Outside this regime we 
find that Z±(x) = exp[=F2i(o; + i0)(x — Xq)]Z±(xq) is 
the solution of Eq.([ll]). Physically it is clear that expo- 
nentially growing solutions are forbidden, which requires 
Z_(0) = Z + (L) = and implies R = 1 for x < and 
x > L. We conclude that at the boundaries 



ipo = 





(14) 



where Z + (0) and Z^(L) are determined by ip(L) = 
S(L)jj . This implies Z+(0) = ~S 12 {L)/ Su(L). Be- 
cause the matrix elements of S(L) depend on the disor- 
der, the initial vector ipo is stochastic. Note that in text- 
book discussions of multiplicative stochastic differential 
equations one often assumes deterministic initial condi- 
tions IQ . After simple algebra we obtain for the second 
component of ipix) — S(x)^>q in the interval < x < L 

R(x) = S 22 (x) - S 21 (x)S 12 (L)/S n (L) . (15) 

In general we have to rely on approximations to calcu- 
late the time-ordered exponential in Eq.(13). However, 
there are two special cases where S(x) can be calculated 
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exactly. The first is obvious: if A(x) = Ao is indepen- 
dent of x, our spin Hamiltonian H(x) is constant, so that 
the time-ordering operator is not necessary. Eq.(|l5|) can 
then be evaluated exactly for arbitrary L [pjj . If we take 
the limit L — > oo holding x/L fixed, we recover the well 
known square root singularity at the band edges, 



A = £/L, and v — Aq/{2Kq^), Numerical results for 
(p(x,Q)) are shown in Fig.[l]. Due to symmetry with re- 
spect to x = L/2 the local DOS assumes an extremum 
at x — L/2, which in the limit L — > oo approaches either 
zero or infinity. Using the fact that at x = L/2 the cosh 
in the numerator of Eq. (|l^) is unity, we obtain 



e(^ 2 -Ag)H 

hm p(x,uj) = -=-y 



A 2 



0<x/L<l. (16) ( p (L/2,0)) = 



There exists another, more interesting limit where S(x) 
can be calculated exactly. Obviously, at w = the direc- 
tion of the magnetic field in our spin Hamiltonian (^Tj) is 
constant. Although in this case H(x) is ^-dependent, we 
may omit the time-ordering operator in Eq .(|13| ). After 
straightforward algebra we obtain from Eq. (|15|) 

R(x) = coshL4(x) - B(x)]/cosh[A(x) + B{x)] , (17) 

where A(x) = f*dx'A(x') and B{x) = J^dx'A(x'). In 
Eq.(|l7j) it is understood that R(x) stands for R(x,iO), so 
that pyx, 0) = 7r _1 .R(x). We have thus succeeded to cal- 
culate the local DOS p(x, u = 0) of the FGM at the Fermi 
energy for a given realization of the disorder. The special 
symmetries of random Dirac fermions at u> = have re- 
cently been used by Shelton and Tsvelik |2^] to calculate 
the statistics of the corresponding wave-functions. 

To calculate the disorder average of Eq.([l7]), we intro- 
duce P(x; a, b) = (S(a — A(x))S(b — B(x))} and write 



(R(x)) = 



da 



cosh(q - b) 

dbP(x;a,b) — — . (18) 

v / cosh(a + fe) 



Assuming that the probability distribution of A(x) is 
Gaussian with average A and covariance K(x,x'), the 
joint distribution P(x; a, b) can be calculated exactly. 
The integration over the difference a — b in Eq.([l8|) can 
then be performed, and we obtain for the average local 
DOS (p(x,0)} = n-^Rix)), 



(p(x,0)) 



/CO 
dsexp [-s 2 /(2cr 2 )] 
- OO 

cosh[(3(x)s + A (2x - L)} 



7rV27T(T 2 



cosh[s -I- A Q L] 



(19) 



Here a 2 = d(L), a{x) = 2[C 1 (x)C 2 (x) - Cl(x)]/a 2 , and 
(3{x) = [d(x)-C 2 (x)]/a 2 , with 



C^x) 



dx' / dx"K(x',x") 



(20) 



C*2 (x) is defined by replacing the range of the integrals in 
Eq.(p0[) by the interval [x,L], and Cs(x) is obtained by 
choosing the interval [0, x] for the x'- and [x, L] for the 
cc"-integration. 

We now specify K{x,x') to be of the form (]l|). Then 
tr 2 , a{x), and f3(x) are easily calculated. It is convenient 
to introduce the dimensionless parameters L = 2Kq£L, 



7TV27TCT 2 J-oo COShfs + VL\ 

where cr 2 = L[l - A(l - e" 1 ^)] and 

/(A) = 1 - A[3 - 4e- 1/(2A) + e- 1/x ] . (22) 

As shown in Fig.|^, /(A) is positive and monotonically 
decreasing. Let us first consider the case v = 0. This 
corresponds to the model discussed by Sadovskii |4 1 with 
real A(x). For L 3> 1 the s-integration in Eq.(2l|) is 

easily done, and we find (p(j,0)) oc L -1 / 2 exp[^/(A)]. 
Keeping in mind that for any finite £ the parameter A = 
£/L vanishes for L — > oo, it is obvious that in this limit 
(p(L/2,0)} is infinite. From Eq.(|l9|) it is easy to show 
numerically that this is also true for \imi,^ oc (p(x 1 0)) in 
the open interval < x/L < 1. We have thus proven 
that for L — > oo and arbitrary £ < oo the average DOS 
of the FGM is infinite at the Fermi energy, in agreement 
with general symmetry arguments . 

For finite v a careful analysis Jl7j] of Eq. (|l]) shows that 
there exists a critical value v c (X) such that for \v\ > v c 
the local DOS (p(L/2,0)) scales to zero in the thermo- 
dynamic limit, and a pseudogap emerges. We obtain 



^ c (A) = [l-A(l-e- 1 A)] 1 / 2 [/(A)] 1 / 2 



(23) 



see FigJ|. For A = Eq.@ yields v c (0) = 1. At the first 
sight this seems to contradict the result of OE ||, who 
found pseudogap behavior already for \u\ > 1/2. One 
should keep in mind, however, that we have set lo = be- 
fore taking the limit L — > oo, while in Ref. Q these limits 
are taken in the opposite order. The non-commutativity 
of these limits is well known from the calculation of the lo- 
cal DOS of the Tomonaga-Luttinger model with a bound- 
ary |21|. Interestingly, for K(x,x') = 2D6(x — x') there 
exists a regime 1/2 < v < 1 where in the thermodynamic 
limit is discontinuous at lu = 0. This follows from 

the fact that according to OE (p(cj)) oc |w| 2iy_1 — * 
for cij — > 0, whereas we have shown that (p(0)) = oo. 
For v — and large but finite £ we conjecture that the 
frequency-dependence of limL^ 00 (p(ti;)) is qualitatively 
similar to the behavior found by Fabrizio and Melin 111]] 
for random Dirac fermions with a special type of disorder 
p2j. Specifically, we expect that for frequencies exceed- 
ing a certain crossover frequency oj* the average DOS of 
the FGM shows pseudogap behavior, which is correctly 
predicted by Sadovskii's algorithm M. However, this al- 
gorithm misses the Dyson singularity, which emerges for 
frequencies \u>\ lu* for any finite value of £. 
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In summary, we have developed a non-perturbative 
method to calculate the Green's function of the FGM. 
Our main result is the proof of the existence of a singular- 
ity in (p(0)} for any finite value of the correlation length £ 
as long as A(x) is real and |Ao| is sufficiently small. For 
Ao = we have shown that with open boundary con- 
ditions (p(f,0)) oc exp[Ko£Lf(£/L)], where /(0) = 1. 
Moreover, if we let L — > oo keeping x/L E (0, 1) fixed, 
(p(x,0)) exhibits a similar singularity Q. Thus, disor- 
der pushes states from high energies to the band center. 
For finite Ao this effect competes with the suppression of 
the DOS due to long-range order. In the incommensurate 
case (where A(x) is complex and in Eq. we should 
replace io^A — » rr + A — rr_A* Jj it is known that (p(0)) is 
finite in the white- noise limit jl2],^3| . We expect that this 
remains true for arbitrary £. In fact, there exists numeri- 
cal evidence that in the incommensurate case (p(^)) can 
be accurately calculated from Sadovskii's algorithm [E3 . 
For a comparison with experiments one should keep in 
mind that any violation of the perfect charge conjugation 
symmetry will wash out the singularity at u> = 0. It is 
therefore unlikely that the singularity is visible in realistic 
materials, although an enhancement might survive. The 
fact that the singular behavior of (p(0)) in the FGM with 
real A(x) is only destroyed if \v\ = | A |/(2iYo£) exceeds a 
finite critical value implies that in commensurate Peierls 
chains true pseudogap behavior should emerge gradually 
below the Peierls transition when the order parameter Ao 
is sufficiently large. Lee, Rice, and Anderson |3j came to 
a similar conclusion within perturbation theory. 
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FIG. 1. Disorder average (R(x,i0)) — n(p(x, 0)) for L = 4 
and A = as function of x/L, see Eq.(|l9|), From top to 
bottom: v = 0, 0.5, 0.75, 1, 1.25, 1.5, 2. 
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FIG. 2. Plot of /(A) and u c {\) denned in Eqs.Q and 
(§)• 
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